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ABSTRACT 

The Laughlin states for N interacting electrons at the plateaus of the fractional Hall 
effect are studied in the thermodynamic limit of large N . It is shown that this limit leads 
to the semiclassical regime for these states, thereby relating their stability to their semi- 
classical nature. The equivalent problem of two-dimensional plasmas is solved analytically, 
to leading order for N — > oo, by the saddle-point approximation - a two-dimensional exten- 
sion of the method used in random matrix models of quantum gravity and gauge theories. 
To leading order, the Laughlin states describe classical droplets of fluids with uniform 
density and sharp boundaries, as expected from the Laughlin "plasma analogy". In this 
limit, the dynamical W^-symmetry of the quantum Hall states expresses the kinematics 
of the area-preserving deformations of incompressible liquid droplets. 
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Introduction 

The current theory for the explanation of the plateaus in the fractional quantum Hall 
effect |1J is based on the seminal work of Laughlin The main idea is the existence of 
incompressible quantum fluids at specific rational values of the electron density. These are 
very stable, macroscopical quantum states with uniform density p(x) = v eB /he = const., 
v = 1/m, m = 1,3, 5, . . ., which possess an energy gap (here B is the external magnetic 
field). Incompressibility accounts for the lack of low-lying conduction modes, which causes 
the longitudinal conductivity a xx to vanish, while the overall rigid motion of the uniform 
droplet of fluid gives the rational values of the Hall conductivity a xy = ue 2 /h. 

In Laughlin's theory, the incompressible quantum fluid configurations of iV interacting 
electrons are described by the wave functions 

1 N 

n ~^pJ2\ z *\ 2 

i/> m (zu . . . , z N ) = C j [ (Zi-Zj) m e i=i , m = 1,3,5,..., (1) 

i<j=i 



where £ = \j2UcjeB is the magnetic length, C = £- Ar -™- Ar ( Ar - 1 )/ 2 ; aric [ {zi,zj} are the 
coordinates of the particles in complex notation. The many-body properties of these states 
are described by expectation values of the form 



1 r N rl 2 

^) = i- /n^ ^iw 2 . ( 2 ) 



where 



Zm. .1 M ^ 2 
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These are interpreted as averages in the reduced statistical problem of a two-dimensional 
one-component plasma [[|, characterized by the effective temperature T oc 1/m. This 
"plasma analogy" has been an important source of physical insight. For example, the 
property p(x) = const, of the electron ground states has been deduced @ from the fact 
that the plasma is a liquid at high effective temperatures, m <^ 70 [Q. Debye screening, 
translational invariance of the liquid and numerical results on correlation functions have 
been used to derive the ground state energy , and the properties of quasi-particle 
excitations (like their fractional statistics and their energy gaps 0). 

A challenging problem is to devise a reliable analytic method for computing the plasma 
partition function (^). Similar expressions for one- dimensional plasmas are found in the 
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theory of random matrices 0. Actually, Laughlin's plasma (^) has been recently refor- 
mulated as a two-dimensional matrix model ||. Given that matrix models can be solved 
with the 1/iV-expansion technique ||, we are lead to conclude that the same technique 
might be successfully applied to Laughlin's theory It is the purpose of this letter to set 
the frame for such a large iV expansion in Laughlin's theory. 

We derive the leading term by the saddle-point approximation of the plasma partition 
function (H). This correctly describes the semiclassical incompressible fluid state of the 
Laughlin plasma analogy. Our main physical point is to show explicitly that, for the 
Laughlin states, the thermodynamic limit N — > oo implies the semiclassical limit % — ► 0, 
and viceversa. While this equivalence between limits is well-known in gauge theories and 
matrix models, it acquires a different physical status in our problem, because N is a not 
a tunable parameter but, rather, it is naturally forced to take large values. Therefore, the 
stability of the Laughlin fluid follows from its prominent semiclassical nature*. 

In the second part of this letter, we discuss the W^- symmetry underlying both Laugh- 
lin's incompressible fluids and the c = 1 matrix models. This symmetry has been explicitly 
shown to account for the incompressibility of the v = 1 Hall ground state [|IT] (similar argu- 
ments apply also to the v = 1/m Laughlin states |[11|| ). Here we show that, in the large N 
limit, the Woo transformations reduce to the classical deformations of liquid droplets which 
preserve the area. These deformations, called area-preserving diffeomorphisms, satisfy the 



Woo-algebra [0**. This classical droplet picture has been already developed for the c = 1 
matrix models of string theory in Refs. [[[4] ||15| , where the relation to the Landau level 
problem has been also recognized. By identifying the H — > and the iV — > oo limits, we 
show that this picture nicely fits into Laughlin's plasma analogy. 



Saddle-point approximation 

Let us start by recalling some known facts about the plasma partition function Z m 
in eq.(|3|). For m = 1, it can be computed exactly by using the orthogonality of the first 
Landau level wave functions, 



1 



>,\ 2 /2t 2 



(4) 



* A similar conclusion was reached in [1C] by a functional approach. 
** See also [13] for a related discussion of classical incompressible fluids. 
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where k is the angular momentum eigenvalue. The result is 
logZi = log | TV! Y[ n\] 

\ n=0 / ( 5 ) 

N 2 ( 3\ N 1 

= _h gJV--J + log AH + -log27r - ^logiV + 0(1) . 

We are interested in the observable one-particle density, 

pi(x) = (fii|*t( x )*(x)|0 1 ) , (6) 
where is the v = 1 ground state and \& the field operator in fermionic Fock space, 

no 

k=0 

The density is easily computed: 



V(z,z) =Y J F kVk(z,z), {F kl F}} = 5 Kl . (7) 



i Ar_1 i ik 



£ 2 7T ^ fc! 

k=0 

and is plotted in Fig. 1 for N = 50. It is constant for r <C £y/N, and drops rapidly to zero 
around r ~ £\fN . This is the density profile of a quantum droplet of incompressible fluid. 
The fluid character is reflected by the uniform value in the interior, and incompressibility 
follows from the gap for density fluctuations (the cyclotron energy). Quantum behaviour 
is apparent from the smooth boundary, where the occupation probability is neither zero 
nor one. In contrast, the density of a classical droplet of liquid has a sharp boundary. 

Actually, the large N limit of the quantum density (Q) is a classical density. Let us 
introduce the rescaled coordinate 

w = ~4= , (9) 



and the rescaled density pi(w) = pi(\z\ = VW\w\) , satisfying J d 2 w pi(w) = 1 . The 
rescaled density possesses a finite large N limit, 

w to P,(N = vSH) = ^e(i-^£) , (io) 

where is the step-function. The sharpness of the boundary is a first indication of the 
classical nature of the iV — ► oo limit. In eq. (|TH|), £ 2 has to be understood as the classical 
parameter setting the scale for the electron density through 

N B B 1 
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for uniform filling v and given area A of the sample. 

Equation (|TT]) implies that the naive semiclassical limit h —* with all the other 
parameters fixed cannot be taken in our problem. For a given external magnetic field B 
and type of fluid characterized by u, the limit % — > enforces N — ► oo if the system is to 
have a finite macroscopical area A. This is another general argument for the advertised 
equivalence of the N —>■ oo and ti — > limits. 

No analytic expression is known for the densities p m , m = 3, 5, . . ., corresponding to 
the Laughlin states (|l|). Numerical studies |16|] [17|] for large number of particles (up to 
N = 200), show a constant density p m = l/mni 2 for r <C iVmN , followed by an 
upward bump near the boundary region, where p m drops rapidly to zero. The fact that 
this bump does not seem to decrease rapidly when the particle number iV is increased up 
to N = 200 has led the authors of jH]] to conclude that Laughlin's wave function does not 
describe a uniform quantum fluid. Actually, we are now going to show that this is not the 
case. Indeed, the limiting large N form of the densities p m will be shown analytically to 
be: 

j limp ro (N = ViVH) = -^0^-1^). (12) 



These describe again classical droplets of incompressible fluid; therefore, the bump seen in 



ref. |T7|] has to disappear for iV — > oo. 

This result can be obtained by extending the saddle-point technique of Brezin, Itzyk- 
son, Parisi and Zuber || to Z m in eq.([$D*. We first rewrite: 



r N d 2 

Z ™ = / II — ^ exp(-NH m [w]) 
J i=i 71 



AT I I o N 
tj r i Fi m \^ i \Wi-Wj\ 

i=l i<j = l 



(13) 



For large N, the particles are driven into a saddle point configuration {wi = w®}, deter- 
mined by the equation 



By considering a lattice decomposition of the plane, we can replace the sum over particles 
with the sum over cells times the characteristic function for cell occupation. For iV — » oo, 



See also [18] for a similar approach. 
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the latter becomes a continuous distribution for the rescaled variable w, which equals the 
electron density p m {w) to leading order. Therefore, we perform the replacement 

-> J d 2 z P m(z) = N J d 2 w p m (w) . (15) 



The saddle-point equation (|TJ) becomes the integral equation 



1 2 I i2 I Pm( w ) , , ... 

w = m £ a w , (16) 

J w — w' 

whose solutions are subjected to the normalization 

J d 2 w p m (w) = 1 . (17) 

The double integration makes the integral equation (|T6|) more involved than the cor- 
responding one for one-dimensional matrix models |§ - no general solution is known to 
us. However, it is easy to check that the rotational invariant density (|12]) is a solution***. 
Furthermore, the following argument shows that this solution is stable, i.e., it is a local 
minimum of H m in eq.(^3|). By scaling the variable w = \frnq, one obtains the identity 

Z m = exp ^^ mlogm) (Nmf J e~™ NH ^ , (i 8) 



which can be approximated semiclassically to second order, yielding 



where Vi = (u>i, Wi) . This equation shows that quadratic fluctuations are independent of 
m. Moreover, for m = 1 the saddle-point solution is clearly quadratically stable, since it 
agrees with the exact solution (^)(^0|). Therefore, the solution ([12]) is stable for any m. 

The value of the Hamiltonian at these saddle points can also be easily computed in 
the continuum approximation and it reads: 

N 2 m ( S\ 
logZ m = -— logiVm - + log AH + ... , (20) 



2 V 2 



** The point w = w' excluded in the sum (|T^) causes no harm in the following integral for 
continuous p functions. 

*** To see this, it is easier to perform first the angular integration using the theorem of residues. 
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which also matches smoothly the exact value for m = 1 in eq.([5]) to leading order. 

The occurrence of other stable solutions of lower "energy" than ([20]) for m = 3, 5, . . . , 
is not excluded by our analysis, but it is very unlikely for small values of m. Numerical 
simulations of the two-dimensional plasma || indicate a phase transition to a Wigner crys- 
tal only at a large value m cr it ~ 70. It would be very interesting to find the semiclassical 
solution with broken rotational invariance corresponding to the Wigner crystal. Let us 
also remark that the one-dimensional plasma of matrix models has only one phase 0. 

In conclusion, for N — > oo we confirm that the Laughlin wave functions ([D describe 
droplets of uniform fluid with densities p m = 1/miil 2 and sharp boundaries. The inter- 



esting structure at the boundary of the droplets found by the numerical calculations [|17 



might be related to edge excitations [19], which are subleading 0(1/N) boundary effects 
20. 



Observables to leading order 

Some physical information on quasi-particle excitations can be obtained by evaluat- 
ing observables (|2]) within the saddle-point approximation. More precisely, all Laughlin's 
results based on the plasma analogy can be rephrased in this approximation. Following 
the review article || we can, e.g., verify the normalization of the wave function for one 
quasi- hole at the point z, 

N 

i>QH(z; zi, . . . , z N ) = e" |z| /w Y[(z - Zi) ip m (zi, . . . , z N ) = S[z, Zi] ip m (zi, . . . , z N ) , 

i=l 

(21) 

with ifj m the Laughlin wave function ([!]). Let us compute H^QhII 2 , i.e. eq.(0) for the 
operator O = S. Repeating the previous steps for the modified plasma with one added 
charge, we find that the saddle-point equation ( |16|) is not modified to leading order, and 
that the saddle-point value of this observable is indeed unity. Actually, the saddle-point 
solution ( |i~2|) is not modified by the inclusion of any finite number of charges. 

The wave function of two quasi-holes can be treated similarly. One finds that the 
correlation among them, showing their fractional statistics, is a subleading effect. This 



agrees with the result of the theory of edge excitations [|19] , in which it has been shown 
that quasi-particle correlations are of order 0(1/N) |f20"[ . As expected, quasi-particles, 
which are quantum effects, are not seen to leading large iV order. 
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Droplet picture and symmetry 

So far, we have been considering the large iV limit of the density from a computational 
point of view. Now, we would like to discuss the geometrical interpretation of this limit. 
Before doing that, let us recall the Woo dynamical quantum symmetry of the v = 1 



ground state recently found in ref. |Tl[. There, we constructed operators £ n ,m, living in 
the first Landau level, 

N 

£ n , m = (4) n+1 K + \ n,m> -1, (22) 
i=i 

where bi, b\ are the harmonic oscillators for angular momentum J excitations, 

hi = § + £di , b\ = | - lh , [h, b]} = , J = b h ■ (23) 

i 

The £ n , m satisfy the Woo algebra* 

( Minim, k) . . . , . \ 

f-< (m - s)\{k - s)\{s + 1)\ J 

= ((to + l)(k + 1) - (n + + 1)) £ n+ fc, m +Z + (■ • ■) C n +k-l,m+l-l + 

(24) 

In particular, the angular momentum is J = Cqq. The subalgebra 

[£oo, £n,m] = (n - m)C njm , (25) 

shows that the C n ^ m are raising (n > m) and lowering (n < m) operators for angular 
momentum. 

The quantum symmetry of the v = 1 ground state is encoded in the invariance of 
the ground state under an infinite (for iV — > oo) set of these transformations, i.e., the 



highest- weight conditions [11 



£n,m i>i(zi, . ■ • ,z N ) =0, n<m. (26) 

These are interpreted as the algebraic conditions of incompressibility, since eq.(|26|) means 
that all transitions lowering the angular momentum of the ground state, i.e., compressions, 
are impossible. 



* More precisely, this is a V^i+oo algebra. 
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Next, we discuss the corresponding classical picture. As shown before, when iV — > oo, 
the quantum density of the Laughlin states at v = 1/m reduces to the profile of a classical 
droplet of liquid with uniform density p m = 1 / rmri 2 . Consider now a deformation of this 
droplet. The density cannot change locally, due to incompressibility. Moreover, its space 
integral gives the particle number N, and is, therefore, constant. Thus, the area occupied 
by the droplet stays constant, i.e., deformations can only produce droplets of the same area 
and different shapes. Therefore, the different configurations of the incompressible fluid are 



related by area-preserving diffeomorphisms, whose generators satisfy the Woo algebra ||12|| . 
In the following, we give an explicit derivation of the action of on the observable 
density and its classical limit Woo, thereby confirming this picture. 

To this end we use a different basis for the generators of Woo- By using this new 
basis, we will stress the analogy to the parallel discussion of the % — > limit of the c = 1 



matrix model in (1 + 1) dimensions, as formulated in [H][I5|, which possesses the same Woo 
symmetry. The basic reason for this analogy is that the Hilbert spaces of the two systems 
are isomorphic. In the matrix model, one considers the Hamiltonian H = \ (p 2 — x 2 ) 
yielding a real representation of the harmonic Fock space. In the first Landau level, we 
have the holomorphic representation (|23|), where angular momentum plays the role of the 
Hamiltonian. As is well known, in Bargmann space the two coordinates (z, z) become 
conjugate variables of a (1 + 1) -dimensional phase space [ pTf . 

In analogy with [[TJJ, we consider the Wigner phase-space distribution: 

W(k,k) = J d 2 z^{z)e ¥ ^ +kb ) = e £2k ~ k / 8 J d 2 z ^ (z) e +^ z+ks) tf(s) . (27) 

This is the generating function for the operators* C n ,m in eq.(|22"|), 



n,m=0 \ / \ / 

where : : denotes Weyl normal ordering |THf , e.g., : (b^) 2 b := | (^(b^) 2 b + b^bb^ +b(b^) 2 ^j. 
This shows that the generating functions W(k, k) form a different basis for the gen- 
erators. In terms of these new generators, the algebra (p4|) acquires the compact 
form 

[W(k,k),W(p,p)} = 2smh(^-(pk-pk)] W(p + k,p + k) . (29) 



* In equation (E7]), \P is the field operator (^); the £ n ,m appearing hereafter are, therefore, 
expressed in the second quantized formulation JO ]. 
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This result is analogous to the (1 + l)-dimensional result of [H], due to the previously 
noticed mapping between Hilbert spaces. 

The Fourier transform of W(k, k) is the one-particle density operator of eq.(|6|), 



p(z,z) 



d 2 k 



e -i(kz+kz) W ( kjk j e 



-£ 2 kk/8 



(30) 



apart from a normal-ordering factor. We recall that, originally Wigner distributions 
were introduced as the quantum analogs of phase-space distributions of classical statistical 
mechanics. This is because quantum expectation values can be written as classical averages 
in phase space with these distributions. In our Landau level problem, the phase space 
is given by (z, z); thus a phase space distribution is actually a two-dimensional space 
distribution, as in (|30|). 

The action of a transformation generated by W(k, k), with infinitesimal parameter 
e(fc, /c), on the ground state density (||) is given by 



S £ pi(z, z) 



d 2 k e(k,k) 5 k t pi(z,z) 



(31) 



where 



8 k - kP i(z,z)= i(Qi\\puW(k f k)]W 



£ 2 kk 



= ze 



U 2 k^- 
e 2 dz 



e 2 8z 



-{kz + kz) 
e 2 pi(z,z) 



(32) 



This expression becomes more transparent in the large iV limit, which is achieved, as before, 
by introducing the rescaled coordinate w = z/yN, and correspondingly, the rescaled 
momentum k = kyN. After rescaling, eq. (|3^) reads, 



5k,r pi(w,w) = 



PuW 



K K 



N VN 



|Oi> 

a 2 _ d \ i 



Pkk I i£_ d_ 

= ie N~8 e 2N K dw - e 2iV* ' dw 



X KW + KW) 

e 2 Pi(w, w) 



J 



where we recall that £ 2 /N = 2hc/eBN. In the large A" limit, eq. (|33]) reduces to 



(33) 



PB 



(34) 



where 6)$ = f S K>h , and {f,g} PB = -| ((df/dw)(dg/dw) - (df/dw)(dg/dw)), with 
B = eB/2c, denotes the correct Poisson bracket with respect to the classical variables 
w and w. This equation shows that, to leading large N order, W {^k/^/N ,k/\/N j is 
the generating function of canonical, and therefore area-preserving, transformations in 
the two-dimensional phase space (w,w). Eq. (|34] ) makes manifest the classical nature 
of a Woo-transformation acting on the density, in the large N limit. As we now show, 
in this limit the quantum algebra reduces to the classical algebra of area-preserving 



diffeomorphisms u^. Our result matches the h — > limit of Refs. [14] [15], thereby further 
confirming that the iV — > oo is a semiclassical limit. 

We now verify the Woo algebra for the classical analogs of the operators W(k, k) 
and £ n} m acting by Poisson brackets on functions of the holomorphic phase space. To 
this end, we perform the coordinate and momentum rescalings as before, and we use 
the correspondence between operators and classical functions, N W ( k/ y/N , R/VW 



) and [ , ] — > ih{ , }pb- The classical limit of eq. (|29|) reads 
|h^(k,k), W (cZ) (A,A)}^ = (Xk- Xk) W (c1 \k + X,R + A) , (35) 

which verifies the identification W^ c1 '(k, R) = ~ e i( KW + K ' w ) from eq. ([34]). This, in turn, 
identifies the classical limit of the operators C n ,m through eq.(|2^), 



/ \ n + m 

\^=) Cn,m - 41 = ~W n+1 W™ +1 . (36) 



Their classical algebra is 

(4% £ ( $}pb = -~ ((m + l)(k + 1) - (n + + 1)) C^ k>m+l , (37) 

which agrees with the classical limit of eq. (|24]) . Equations ( |55j ) and ( |3"TD are equivalent 
forms for the it^ algebra of area-preserving diffeomorphisms. Note that the classical 
functions Cn%, and and their algebras ( |35D and (l37| ) can be also derived directly 



from the canonical treatment of the classical theory describing the dynamics of the first 
Landau level, the "topological quantum mechanics" of refs. |2^] (see also [|TT|). 



Having identified the classical Cnfm, we can evaluate their action on the classical 



density ( |1Q|) by using eq. (|3 



6W Pl (w,w) = { p x {w,w) , w n+1 w™ +1 } pB 

pn+m-2 / |,„|2\ (38) 

= i{n _ m ) — -— e^-^ d 5 ( 1 - N 



nB \ £ 2 
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where 9 is the angular variable on the circle delimiting the classical droplet. These varia- 
tions correspond to density waves localized on the one- dimensional sharp boundary of the 
classical droplet. The quantization of these edge waves leads to a (1 + 1) -dimensional c = 1 
conformal field theory |TI| |20| . 

So far, our explicit derivation of the classical Woo algebra referred only to the 
v = 1 ground state; however, it can be easily extended to the Laughlin states at 
v = 1/m. Indeed, their classical densities p m are related to p\ by the simple scal- 
ing p m [yS) = (l/m)pi(w/y/m). Therefore, we can apply the same scaling in the classical 
phase space to obtain the corresponding action fl34D -(|37p for v = 1/m. We can also 
argue that the dynamical quantum symmetry of the Laughlin states must be the quan- 
tum extension Woo. Actually, the analogs of the symmetry conditions (^6|) were found in 



Ref. [Tl|]. However, in going from to W^, one needs the explicit form of the higher 
order corrections 0(1/N k ), k > in the algebra W^, as in eq.(|24|). These are not unique 
and their correct form is not presently known to us for v = 1/m. 
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Figure caption 
Fig. 1 

The density profile in units of l/n£ 2 for the first Landau level filled up to L = 50 as 
a function of r/t. 
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